Casimir-Lifshitz force out of thermal equilibrium 
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We study the Casimir-Lifshitz interaction out of thermal equilibrium, when the interacting ob- 
jects are at different temperatures. The analysis is focused on the surface-surface, surface-rarefied 
body, and surface-atom configurations. A systematic investigation of the contributions to the force 
coming from the propagating and evanescent components of the electromagnetic radiation is per- 
formed. The large distance behaviors of such interactions is discussed, and both analytical and 
numerical results are compared with the equilibrium ones. A detailed analysis of the crossing be- 
tween the surface-surface and the surface-rarefied body, and finally the surface-atom force is shown, 
and a complete derivation and discussion of the recently predicted non-additivity effects and new 
asymptotic behaviors is presented. 
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I. INTRODUCTION 

The Casimir-Lifshitz force is a dispersion interaction of 
electromagnetic origin acting between neutral dispersive 
bodies without permanent polarizations. The original 
Casimir intuition about the presence of such a force be- 
tween two parallel ideal mirrors fTl (or between an atom 
and a mirror, i.e. the so called Casimir-Polder force H) 
was readily extended to real materials by Lifshitz [1, 0,111 ■ 
He used the theory of electromagnetic fluctuations devel- 
oped by Rytov [6| to formulate the most general theory of 
the dispersion interaction in the framework of the statis- 
tical physics and macroscopic electrodynamics (see also 
0). The Lifshitz theory is still the most advanced one; 
today it is extensively accepted providing a common tool 
to deal with dispersive forces in different fields of science 
(physics, biology, chemistry) and technology. 

It is useful to stress here that the geometry of the sys- 
tem is relevant for the explicit calculation of the force, 
but does not affect the nature of the interaction that pre- 
serves all its peculiar characteristics and relevant length- 
scales. For this reason we refer to the Casimir-Lifshitz 
force for all geometrical configurations. In particular, in 
this paper we are interested in the force between flat and 
parallel surfaces of two macroscopic bodies, and between 
a surface and an individual atom. 

The Lifshitz theory is formulated for systems at ther- 
mal equilibrium. In this theory the pure quantum effect 
at T = is clearly separated from the finite temperature 
effect. The former gives a dominant contribution at small 
separation (< 1 /im at room temperature) between the 
bodies and was readily confirmed experimentally with 
good accuracy [see Q (surface-atom), [1, 

na, El [11 

(surface-sphere), [ll] (surface-surface)]. 



The thermal component prevails at larger distances 
and was measured only recently at JILA in experiments 
with cold atoms [l^. These experiments are based on 
the measurement of the shift of the collective oscillations 
of a Bose-Einstein condensate (BEC) of trapped atoms 
close to a surface The JILA group measured the 

Casimir-Lifshitz force at very large distances (~ fO^m) 
and for the first time showed the thermal effects of the 
Casimir-Lifhitz interaction (and indeed of any dispersion 
interaction), in agreement with the theoretical predic- 
tions • This measurement was done out of thermal 
equilibrium ^18^, where thermal effects are stronger. 

There was an interest in configurations out of ther- 
mal equilibrium since the work by Rosenkrans et al. 
[1^ (atom- atom). Surface- atom interaction was ana- 
lyzed by Henkel et al. (20| and by Antezza et al. 
(itI [21I [2^ m, [23] . Surface-surface force was investi- 
gated by Dorofeyev et al. [1^ HI] and Antezza et al. 

23 , [13] • For a review of non-equilibrium effects see also 

23. 



'Electronic address: ' antezza@science . unitn. itI 



Further non-equilibrium effects were explored by 
Polder and Van Hove [1^, who calculated the heat-flux 
between two parallel plates, and Bimonte [2^, who ex- 
pressed fluctuations of fields for the metal-metal config- 
uration in terms of surface impedance. 

The principal interest in the study of systems out of 
thermal equilibrium is connected to the possibilit y of tun- 
ing the interaction in both strength and sign [ItI H^. 
Such systems give also the way to explore the role of ther- 
mal fluctuations, usually masked at thermal equilibrium 
by the T = component which dominates the interaction 
up to very large distances, where the actual total force 
results to be very small. 

A crucial role in explaining the peculiarity of the non- 
equilibrium surface-atom force is played by cancellation 
effects between the fluctuations of the different compo- 
nents of the radiations, as the incident to and emitted by 
the surface [I7i ]. 
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In this paper we present a detailed study of the 
Casimir-Lifshitz force out of thermal equilibrium, with 
particular attention devoted to the surface-surface and 
surface-atom interactions. We perform a systematic in- 
vestigation of the contributions to the force coming from 
the propagating and evanescent components of the elec- 
tromagnetic radiation. The large distance behaviors of 
these interactions are extensively discussed, both analyt- 
ically and numerically, and comparisons with the equi- 
librium results are done. We perform a detailed analy- 
sis of the relation between the surface-surface interaction 
when one body is rarefied (surface-rarefied body force) 
and the surface-atom force. We also present a complete 
derivation and discussion of the recently predicted non- 
additivity effects and new asymptotic behaviors noted in 
[23,]. 

We are interested in the force occurring between two 
planar bodies, which are kept at different temperatures 
and separated by a distance I. We consider that the bod- 
ies are thick enough, in order to exclude possible effects of 
the presence of the vacuum gap on the radiation outside 
the two bodies. We also assume that each body is in local 
thermal equilibrium, the whole system being in a station- 
ary state. In our configuration the left-side body, 1, has 
a complex dielectric function ei(ijj) — e[{uj) + ie'^uj), 
occupies the volume Vi and is held at temperature Ti. 
The right-side body, 2, has a complex dielectric function 
£2(w) = £2('^) + *£2('^)j occupies the volume V2 and is 
held at temperature T2. First we assume that each body 
fills an infinite half-space, in particular Vi and V2 coincide 
with the left and right half-spaces, respectively. Later 
we consider a more general situation of two parallel thick 
slabs with the external regions shined by the thermal 
radiations at arbitrary temperatures. In this case addi- 
tional distance- independent contributions to the pressure 
are present. Finally we will consider the case in which 
one of the two bodies is rarefied. In this case the in- 
terplay between the finite thickness of the body and the 
non equilibrium configuration leads to different interest- 
ing behaviors of the pressure. 

The general problem can be set in the following way, for 
two bodies occupying the two half-spaces. Let us choose 
the origin of the coordinate system at the boundary of the 
half-space 1 and let us set the z-axis in the direction of 
the half-space 2 [see FigH] . The electromagnetic pressure 
between the two bodies along z can be calculated as [lO, 

Hi 

P"'='i(Ti,T2,0 = (r,,(r,i)), (1) 

that should be regularized by subtracting the same ex- 
pression at separation I 00. In Eq.(IT]), r is a generic 
point between the two bodies, and 

r,,(r,i) = - ^ [E^{r,t)Ef3ir,t) + B^{r,t)Bfj{r,t)] , 

(2) 

is the zz component of the Maxwell stress tensor in the 
vacuum gap. Here Aq,^ is a diagonal matrix with An = 
A22 = 1 and A33 = -1. 




FIG. 1: Schematic figure of the surface-surface system out 
of thermal equilibrium. Here the two bodies occupy infinite 
half-spaces. 



To calculate the pressure (P) one must average over the 
state of the electromagnetic field the squares of the spa- 
tial components of the electric and magnetic field E(r, t) 
and B(r, t), which appear in Eq. 

Before starting with the analysis of the problem we 
mention the structure of this work in the following out- 
line. In Sec. |TT]we develop the formalism, introduce the 
role and the description of the fluctuations of the electro- 
magnetic field and specify the approach we adopt to deal 
with the surface optics. In Sec. IIIII we recall the main 
results of the surface-surface Casimir-Lifshitz interaction 
at thermal equilibrium, and in particular specify the dis- 
tinction between the T — (purely quantum) and the 
thermal contribution to the force, generated by the ra- 
diation pressure of the thermal radiation. In Sec. IIVI 
we present a detailed derivation of the surface-surface 
pressure out of thermal equilibrium P'^°'^{Ti,T2,l). In 
Sec. |V]we show an alternative and useful expression for 
P™'^{Ti, T2,l), together with numerical results relative to 
particular couples of dielectric materials (i.e. fused silica- 
silicon and sapphire- fused silica). In Sec. IVII we deal with 
the distance-independent terms in the pressure due to the 
finite thickness of the two bodies, and the eventual effect 
of external radiation at different temperature impinging 
the external surfaces. In Sec. I VIII we derive the large dis- 
tance behavior of the surface-surface pressure out of ther- 
mal equilibrium, and discuss the role of the propagating 
waves (PW) and evanescent waves (EW) contributions. 
We also make a comparison with the corresponding terms 
of the pressure at thermal equilibrium. In Sec. IVIIII we 
consider the interaction between a surface and a rarefied 
body and derive the large distance behaviors of the PW 
and EW components. In the same section we stress the 
presence of non additivity in the interaction out of equi- 
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librium (in contrast with the equihbrium case) and show 
the analysis of the crossing between different asymptotic 
behaviors. In Sec. IIXI we show the transition from the 
surface-rarefied body to surface-atom interactions out of 
thermal equilibrium, and demonstrate the essential role 
of finite thickness of the rarefied body. Finally in Sec. |X] 
we provide our conclusions. 

In appendix [X] we give some details on the expression 
of the Green functions we used in our calculation, and in 
appendix IB] we discuss in detail the force acting between 
a surface and a rarefied body of finite thickness. 



II. THE FORMALISM 

Our approach is based on the theory of the fluctuating 
electromagnetic (EM) field developed by Rytov In 
this approach it is assumed that the field is driven by 
randomly fluctuating current density or, alternatively, by 
randomly fluctuating polarization field. In this respect 
the Maxwell equations become of Langevin-type. For a 
monochromatic field in a non-homogeneous, linear, and 
nonmagnetic medium with the dielectric function e(w, r) 
the Maxwell equations become: 

VAE[tj;r]-iA:B[tj;r] = 0, (3) 
V /\B[uj;r]+ike{uj;r)E[uj;r] = -in ik P[uj;r], (4) 

where k — w/c is the vacuum wavenumber, and A is 
the vector product symbol. The source of the electro- 
magnetic fluctuations is described by the electric polar- 
ization P[(jj;r], related to the electric current density as 
J[w;r] = —iujP[uj;r]. We use the following notations for 
the frequency Fourier transforms r] of the quantity 
A(t,r): 

^(i,r)= r°°^e— *4c.;r]. (5) 

To find the solution of the Maxwell equations we use 
the Green functions formalism. A Green function is a so- 
lution of the wave equation for a point source in presence 
of surrounding matter. When this solution is known one 
can construct the solution due to a general source us- 
ing the principle of linear superposition. This method 
takes into account the effects of non-additivity, which 
originates from the fact that the interaction between two 
fluctuating dipoles is influenced by the presence of a third 
dipole. Employing this formalism we can express the 
electric field at the observation point r as the convolu- 
tion 

E[w;r]= jG[uj;r,r'] • P [tj; r'] dr'. (6) 

Here P [uj; r'] is the random polarization at the source 
point r', and G [a';r,r'] is the dyadic Green function of 
the system. Then it is clear that the Green function plays 



the role of the response function in a linear-response the- 
ory. The Green function is the solution of the following 
equation [s^ 

{V A V A -k^e{Lu, r)} G[uj; r, r'] = 4TTk^I6{r - r'), (7) 

where I is the identity dyad. This equation, resulting 
from the Maxwell equations ([3]) and ([4]) and convolution 
has to be solved with proper boundary conditions 
characterizing the fields components at the interfaces, as 
well as the condition required by a retarded Green's func- 
tion i.e. G[uj;r,r'] -> as |r - r'| ^ cx). 

Finally it is useful to recall the relations Gap[uj; r, r'] = 
G/3a[a;; r', r] and G'*^[ijj; r, r'] — Gaf3[—uj;r,r'], that are 
the consequence of the microscopic reversibility in the 
linear-response theory and the reality of the time depen- 
dent fields, respectively. 

A. Field correlation functions 

From Eq.([T]) it is evident that we are interested in the 
time correlations between different components of the 
electric (magnetic) field at equal times. In the quantum 
theory such correlations are described by the averages of 
symmetrized products of the field components: 

(i?„(r,t)i?^(r',i))^y„^ = 

i {E^ir, t)Ep{r', t) + Ep{r' ,t)E^{r, t)) . (8) 

Notice that, although in this paper we are using sym- 
metrized correlations, other possible forms of the corre- 
lation functions could be more appropriate in other situ- 
ations 33] . The correlations ([5]) in terms of their Fourier 
transforms can be presented as 

{E^{v,t)Ep{v\t))^^^^^ 

Using Eq.([6]) these correlations can be expressed via the 
correlations of the polarization field P, which obeys the 
fluctuation-dissipation theorem jsoj ] 

lPc.[Lo-v]Pl[Lo'-r']) = 

\ / sym 

^^-y^ -th (^) 5{u .')5{v r')-5.„ (10) 

expressed via the Fourier transformed P [ti-';r]. Due to 
the presence of the 5 {r — r') factor these fluctuations are 
local. Fluctuations of the sources in different points of 
the material are non-coherent. This permits to assume 
that in the non-equilibrium situation, when temperature 
T is different in different points, the sources correlations 
are given by the same equations. We must emphasize 
that this assumption, even being quite reasonable, is still 
a hypothesis, which is worth both of further theoretical 
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investigation and experimental verification. The prob- 
lem was discussed previously (see particularly [4l|), but 
in our opinion the conditions of applicability of the the- 
ory has not been still established. The same assumption 
was used by Polder and Van Hove [28[ to calculate the 
radiative heat transfer between two bodies with different 
temperatures. 

The assumption (jlOp (local source hypothesis) repre- 
sents the starting point of our analysis allowing for an 



explicit calculation of the electromagnetic field also if the 
system is not in global thermal equilibrium. 

It is now evident that EM field in the vacuum gap is 
given by the sum of the fields produced by the fluctuating 
polarizations in the materials filling respectively the half- 
space 1, with the dielectric function ei{uj) and tempera- 
ture Ti, and the half-space 2 with the dielectric function 
e2(w) and temperature T2. Then the Fourier transform 
of the electric field correlations can be presented as 



sym 



coth 



2M 



coth 



2kBT2 



S {lU — Lo') , 

(11) 



where 5^*^ [i = 1,2) is defined as convolution of two 
Green functions 

5g[c.;r,r']= / dr"G^,[u;;r,r"]G;i,[u;;r" ,r']. (12) 

Here Vi and V2 are the volumes occupied by the left and 
right body, respectively, and the two terms in Eq. (jll[) 
correspond to the parts of the pressure generated by the 
sources in each body separately. 

It is interesting to see how the global equilibrium is 
restored when Ti ^ T2 ~ T in Eq. pTjl . In this case 
Eq. pT|) can be written as 



£;„[w;r]4K;r'] 



^ S {lu — oj') X 
[ drV'(a;,r")G„^[c^;r,r"]G;^[c^;r",r']. (13) 



~ — coth 

2 \2kBT 



The integral over the product of two Green functions is 
connected with the imaginary part of the single Green 
function by the important [4ll |42| relation 

/ dv"e'\ujy)G^,[uj-vy]G;p[oj-v\v'] = 

4rfmG„;3[w;r,r'], (14) 

where f2 is a volume restricted by a surface where the 
Green function vanishes. Keeping in mind that in the 
vacuum gap e" — 0, one can extend the integration 
in Eq. (fT3|) over the the whole space and using (fT4|) 
one recovers the well-known form of the electric fields 
fluctuation-dissipation theorem [3^ valid at a global 
thermal equilibrium: 



I sym 

2^ficoth \MiGo,ti\uj-vy\S{uj-uj') . (15) 



Notice that all fluctuations presented in this section in- 
clude both the vacuum (T = 0) and the thermal fluctu- 
ations. These can be identified with the first and second 
terms, respectively, of the r.h.s. of the identity 



coth 



sign(cj) 1 + 



gRIwI/fcflT _ I 



(16) 



B. The pressure in terms of fluctuations 

The pressure ([T]) can be presented in terms of the 
Fourier transformed fields correlations: 



P'-<i(Ti,T2,0 



1 

8^ 



dcij duj' 



A 



Q/3 



E^[u-v]El[iv'-v'] 



B^[w-v]bI[lo'-t'] 



(17) 



Here the electric and magnetic contributions to the total 
pressure are explicit, and r is a point inside of the vacuum 
gap. The stress tensor is in fact constant in the vacuum 
gap due to the momentum conservation required by a 
stationary configuration (see discussion in section flV A[) . 
In this equation we omitted the symmetrization index 
since the average is taken at the same point r = r'. Using 
Eq.Q it is useful to rewrite expression (|17p in terms of 
the electric fields only [s^l as 



p"'=^(ri,T2,0 



1 

"8^ 



-i{u—uj )t ^ 



2tt 2tt 

Qs. [{Es[u;-r]El[u;';r'])] 



(18) 



Here the pressure is expressed in terms of the correlations 
([TT|) . and the operator 



p-ea^ie^^^a-yC?', j (19) 
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selects the electric and magnetic contributions, given by 
the first and the second term in Eg. p^ . respectively. 

From Eq. (fTHll it is possible to express the total pressure 
as the sum 

P"°^(Ti, T2, 1) = Po(0 + P^:\Ti,T2, 1), (20) 

where the contribution of the zero-point {T = 0) fluctua- 
tions, Pq{1), is separated from that produced by the ther- 
mal fluctuations, P^'^{Ti,T2,l). Furthermore, thanks to 
equation (fTT|) it is possible to express the thermal com- 
ponent of the pressure acting between the bodies as the 
sum of two terms 

pnTi,T2, 1) - p^rm- 0, o + prio, t2,i). (21) 

The pressure at thermal equilibrium P'^'^{T,l), being a 
particular case of ((20| . can be written as 

p°q(T,o-Po(0 + nT(r,0- (22) 

The pressures Po(/) and P°^{T,l) are given by Eq. p^ . 
where the field fluctuations are provided by Eq. after 
the substitution, respectively, of: 

coth ( ) sign(w), (23) 



f hw \ 2 sign(a;) , , 

coth -> ^, V . 24 

\2kBTJ gh\u\/kBT _i V ) 

If one simply perform such substitutions, it is well know 
that equation (|18|) diverges at T = 0, and contains con- 
stant (^-independent) terms in the thermal part. The 
divergence has the same origin as the usual divergence of 
the zero-point fields energy in quantum electrodynamics, 
while the constant terms are related to the fact that we 
consider infinite bodies, and hence we neglect the pres- 
sure of the radiation exerted on the remote, external sur- 
faces of the two bodies. To recover the exact finite value 
for the pressures Po(/), and exclude the constant terms 
in P^^{T, I), one should regularize the Green function in 
the r.h.s of Eq. fTS)) by subtracting the bulk part G^-^-, 
corresponding to a field produced by a point-like dipole 
in an homogeneous and infinite dielectric 0, [H, S^l . In 
fact, the Green function with both the observation point 
r and the source point r' in the vacuum gap (see section 
lA 1[) is given by the sum 

G,, [u;; r, r'] = G^- [u;:, r, r'] + G^/' [lu; r, r'] , (25) 

of a scattered and a bulk term. The subtraction of the 
bulk term corresponds to the subtraction of the pressure 
at Z ^ oo, as prescribed after Eq.(IT|). The expressions for 
the pressure at thermal equilibrium are given explicitly 
in section Hill 

Concerning the thermal pressure out of thermal equi- 
librium P^^'^{Ti,T2j) of Eq. (HOI), it can be obtained 



from Eq.lH]) by using (ITT|) and the substitution ((24)) . 
Also in this case the thermal pressure P^^'^{Ti,T2, 1) con- 
tains an Z-dependent and a constant term, as it hap- 
pens for P^^{T,l) before being regularized. Differently 
from the equilibrium case, here the origin of the constant 
terms is not only due to the absence of the pressure act- 
ing on the remote surfaces, but is also related to the 
fact that out of thermal equilibrium there is a net mo- 
mentum transfer between the bodies. In this case the 
constant terms can remain also after considering bodies 
of finite thickness, and can even be different for the two 
bodies, depending on the external radiations. In the sec- 
tions [Wl and |V] we wiU calculate Pl'^'^iTi,T2,l) for two 
bodies filling two infinite half-spaces, and we will mainly 
discuss the pure /-dependent component. The constant 
terms will be discussed in section IVTl for the general case 
of bodies of finite thickness, with impinging the external 
radiations at different temperatures. 

C. Electromagnetic waves in surface optics 

In this work we formulate the electromagnetic prob- 
lem in terms of s- and p-polarized vector waves and in 
terms of the Fresnel coefficients for the interfaces (37| . 
Such notations are very useful in surface optics. We will 
also employ the angular spectrum representation for the 
description of the EM and polarization vectors. 

If X, y and z are the coordinate unit vectors (with 
real norm equal to 1), one can write the position vec- 
tor as r = R + z z, where the capital letter refers 
to vectors parallel to the interface [R = {Rx, Ry,0)]- 
Let us write the electromagnetic (complex) wave vector 
in the medium m with the complex dielectric function 
e™(w) = £'„(a;) + ie'^{uj) as 

q('")(±) = Q±qi") z. (26) 

Here the sign (-I-) corresponds to an upward-propagating 
(or evanescent) wave, and the sign (— ) corresponds to a 
downward-propagating (or evanescent) wave. The vector 
Q = {Qx,Qy, 0) is the projection (always real) of q('")(±) 
on the interface and the z— component of the wave vector, 
and 

gi") = Vsm k^-Q^ (27) 

is a complex number with a positive imaginary part, and 
positive real part in case Imqi™-* = 0. Real and imaginary 
parts of gi™^ are expressed by the following relations: 

Reqi") = yi[|e„(c.)fc2_g2| + (e;^(^)fc2„Q2)]^ (23) 

Img(™) ^ y i [|e„(c.)fc2 _ Q2| _ (ej^(,,)fc2 _ Q^)].{29) 

Then, if the medium m is non-absorbing (e"^ = 0), 

for Q < e'^ k the wavevector qi™' is real and corre- 
sponds to a wave propagating in the medium m, while 
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for Q > ^y e'^ k the wavevector qi^^ is imaginary and 
corresponds to evanescent wave in the medium m. The 
foUowing identities will be useful: 

2Im<z(™)Reg(™) = e'M. (30) 



It is worth noticing that the wave vectors q'^™)(±) lie in 
the plane of incidence spanned by Q and z. Than one can 
introduce the s- and p-unit complex polarization vectors 



e(™)(±)-QAz, 



(±) =e(™)(±) Aq("'^(±) 



(33) 
, (34) 



that are vectors transversal and longitudinal to that 
plane, respectively. Usually the polarization vector 
ei™^(±) [ep™^(±)] is called transverse electric (TE) 
[transverse magnetic (TM)] since it corresponds to the 
electric [magnetic] field transverse to the plane of inci- 
dence. 

Our geometry consists of two half-spaces labeled with 
m = 1, 2 separated by a vacuum gap. Inside of the gap 
the wave vector q and the polarization vectors e^(±) are 
not labeled and are obtained, respectively, from the def- 
initions ([26]) . (|27| and ([33|) . (pi]) by omitting the apices 
and setting = 1. 

Finally we can introduce the well known reflection 
and transmission Fresnel coefficients for the vacuum gap- 
dielectric interfaces, which for the s- and p-wave compo- 
nents are: 



2gi"'^ 



qzSn 



(35) 



_ 2 ^eraju:) gi""^ 
qi + qz£m{^) 



In particular, the coefficients relate the radiation in 
the vacuum gap impinging the interface m and its part 
reflected back into the vacuum gap. The coefficients tm 
relate the radiation impinging the interface m from the 
interior of the dielectric m and its part transmitted into 
the vacuum gap (see Appendix E]) . 



III. PRESSURE AT THERMAL EQUILIBRIUM 

In this section we briefly recall the main results of 
the pressure in a system at thermal equilibrium. We 
present the thermal component of the pressure as the 
sum of PW and EW components, and in terms of real 
frequencies, which will result useful for the rest of the 
discussion. The results we show for the pressure at 
equilibrium are regularized [see discussion after Eq. 



The Lifshitz surface-surface pressure at thermal equi- 
librium can be expressed in terms of real frequencies as 



xRe 



/ huj 
dQ Q qz g{Q,t^] 



(37) 



where 



'l '2^ 



E [{r^tr'^r^e-^-'^'-lY' ■ (38) 



In the previous equation the multiple reflections are de- 
scribed by the factor 



D 



(39) 



and the reflection Fresnel coefficients for the vacuum- 
dielectric interfaces are defined in Eg. ([55)1 . 

By performing the Lifshitz rotation on the complex 
plane it is possible to write Eq. (|37p in terms of imaginary 
frequencies: 



P°'i(T,0 



kB T 
ksT 



Ax . 



(£io + l)(£20 + l) 



(eio - l)(e2o ■ 
dpp^ gip, i£,n), 



1) 



1 



(40) 



where p = + <?Q^ IQi- The dielectric functions that 
enter to g(j>, i^n) must be evaluated at imaginary frequen- 
cies £1^2 = £i,2(*Cn)i where f„ — 2'KkBTn/h. In the first 
term of we have also introduced the static values of 
the dielectric functions £io ~ £i(0) and £20 — £2(0). 

The pressure at thermal equilibrium includes contribu- 
tions from zero-point fluctuations Po{l) and from thermal 
fluctuations P^^{T, I) as Eq. shows. Po{l) can be ex- 
tracted from I^T]) with the substitutes ([23|) or from ((40|) 
as the limit of continuous imaginary frequency. The final 
result for the T = pressure is 



Po{l) 



27r2c3 



/ d^ / dpp'eg{p,<)- 

Jo J I 



(41) 



The pressure Po(0 admits two important limits, i.e. the 
Van der Waals-London and the Casimir-Polder behav- 
iors, valid at small and large distances, respectively, in 
respect to the characteristic length scale Aopt fixed by 
the absorption spectrum of the bodies (typically is of the 
order of fraction of microns). 

The behavior of the thermal component P^^{T^ I) is 
related to a second length scale, i.e. the thermal wave- 
length 



A, 



he 



(42) 
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which at room temperature is « 7.6 /im. 

Then, the zero-point fluctuations dominate over the 
thermal contribution at small distances Z <C At- In this 
limit behavior of the pressure is determined by the char- 
acteristic length scale Xgpt ^ At- In the interval Xopt <C 
I <C At one enters the Casimir-Polder regime where the 
pressure decays like For distances I ^ Xopt the force 
instead exhibits the l/l^ van der Waals-London depen- 
dence- The possibility of identifying the Casimir-Polder 
regime depends crucially on the value of the tempera- 
ture- The temperature should be in fact sufficiently low 



in order to guarantee the condition At ^ Xopt- 



The last part of this section focuses on the thermal 
component of the pressure, that will be often used along 
the rest of the paper- The pressure P°^{T, I) can be ob- 
tained from (j37p by using (|24P - Since such a component 
of the pressure will be compared with that out of ther- 
mal equilibrium, we show here explicitly its expression 
for PW and EW contributions: 



JO Jk M=s,p ' 



In particular at high temperatures, or equivalently at 
large distances defined by the condition 



I > Xj 



(45) 



the leading contribution to the pressure is given by the 
expression for the total force [3] 



167r/3 



Ax X 



£w + 1 £20 + 1 



(46) 

It corresponds to the first term in Eq- (PD)) and is en- 
tirely due to the thermal fluctuations of the EM field. 
In Ref. [2] the asymptotic behavior (|46|) has been found 
after the contour rotation in the complex w-plane of the 
EW term (|44p , that is partially canceled by the P W term 

One can note that in this regime only the static value 
of the dielectric functions is relevant. The pressure (HS)) 
is proportional to the temperature, is independent from 
the Planck constant as well as from the velocity of light. 
We will call this equation the Lifshitz limit. The pressure 
can be obtained from the thermal free energy T = 
£ — TS of the electromagnetic field (per unit area) ac- 
cording to the thermodynamic identity P = — {dT /dl)j,^ 
where £ and iS are the thermal energy and entropy, re- 
spectively. It is interesting to note that, differently from 
the free energy, the thermal energy £ decreases exponen- 
tially with I, which means that the pressure (|46p has pure 
entropic origin [ssj . 

It is important that at large separations only the p- 
polarization contributes to the force (see for, example in 
[24| . the detailed derivation of the PW and EW com- 
ponents). The reason is that for low-frequencies the s- 
polarized field is nearly pure magnetic, but the magnetic 
field penetrates freely into a non- magnetic material [isj . 



In the limit £10, £20 
two metals 



00 we find the force between 



P™*(T,/) 



8ttP 



C(3). 



(47) 



Let us empathize that this result was obtained for in- 
teraction between real metals [13]. For "ideal mirrors" 
considered by Casimir, both polarizations of electromag- 
netic fields are refiected. In this case there will be an 
additional factor 2 in Eq. (H^ due to the contribution of 
the s— polarization. This ideal case can be realized using 
superconducting mirrors. 

It is useful to note that the surface-surface pressure 
P°'^{T, I) given by the Lifshitz result hides a non triv- 
ial cancellation between the components of the pressure 
related to real and imaginary values of the EM wavevec- 
tors, leading, respectively, to the pro pag ating (PW) and 
evanescent (EW) wave contributions jlTl . [ist . This study 
deserves careful investigation since for a configuration out 
of thermal equilibrium such cancellations are no longer 
present, and the PW and EW contribution will provide 
different asymptotic behaviors. The new effect, as we will 
see, is particularly important if one of the two bodies is 
a rarefied gas. 



IV. PRESSURE OUT OF THERMAL 
EQUILIBRIUM BETWEEN TWO INFINITE 
DIELECTRIC HALF-SPACES 



As was discussed above [see Eqs. ([TT|) and (|2ip ] each 
body contributes separately to the thermal pressure. In 
particular the pressure resulting from the thermal flue- 
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tuations in the body 1 is 



167r3 



Re 



(48) 



where ri = r2 is a point in the vacuum gap and the 
function 5 is defined in Eg. p^ . In Eq. (H51) we used the 
parity properties e"(w) = — e"(— and Ssu [w;ri,r2] = 
Sgj^ [— w; ri, to restrict the range of integration to the 
positive frequencies. It is evident that P^^'^{0, T, I) can be 
expressed similarly to Ea.(|i51). but with £"(ijj) £2('-^) 



and 5* 



1) 



(1) 



Below we specify the expressions of the tensors Sg^ 

and sfj (section llV Ap . calculate the electric and mag- 
netic contributions to the pressure (section lIVBp . and 
finally provide the result for the total pressure in terms 
of P W and EW components (section IIVCP . The total 
pressure will be rewritten in a different form in section 
|V]by using a powerful expansion in multiple reflections. 
In the present and in the next section |V] the pressure is 
calculated for two infinite bodies [see discussion at the 
end of section [TTB] . 



Fourier transforms SgJ [oj; Q, zi, Z2] and s^J [w; Q, zi, Z2] 
one has 



(49) 



By choosing the x-axis parallel to the vector D = Ri— R2 
and defining (j) as the angle between Q and D one gets 
that Qx = Q cos 4>, Qy = Q sin cj) and the polarization 
vectors become 



eW(±) 
e(")(±) 



Here it is evident that |ei'"^(±)p = 1 and 



(I sin 01, — cos(f) sin^/l sin^j, 0) , (50) 
(Tgi™^ cos 0, Tgi"^ sin </., q) . (51) 



eM(±)|^ 



\Sm\ ' 



(52) 



The 5* functions 



In this subsection we show the result for the tensors 
S^gl} and S^^J defined by Eq. In terms of the lateral 

I 



After explicit calculation using the Green function 
given in appendix[X]we find for the s^]^^ and s^^J functions 
the explicit expressions 



s''s'J['^;Q,zi,z2] 



^i(-) ki^¥ ,h l^-l' ' 
^A-KA+) e-^(«^^^+«>^) e2^«^' < + eM,5(-)<,.(-) e-^('^^i-'^>^) e-^i-?^' jr^^^' 



(53) 



E 



^2? 



|e(^)(-)r X 



I 



(54) 



where is defined in Eq. ((39l) . 

It is worth noticing that in the non equilibrium but sta- 
tionary regime the fields correlation functions s'^'^^ are 
not uniform in the vacuum cavity, while on the contrary 



the Maxwell stress tensor T^z (which is related to the 
momentum flux) has the same value in each point of the 
vacuum gap. This is valid also at equilibrium, and is a di- 
rect consequence of the momentum conservation required 
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by a stationary configuration. To show this property one 
can set Zi = Z2 — z in Eq. ([55]) , where the dependence on 
z appears only in the exponential factors [the same would 
happen for Eq. ([5^ ] . Let us note that now the first and 
the last terms in such an expression are proportional to 
g22iing^ g^j^^ g-2zimgj ^ respectively, while the second and 
the third terms are proportional to e"^^^""^^ and 6^^^°'^ ^ 
respectively. As it will be clear in the next section HV Bl 
the first and the last terms will be responsible for the 
PW contribution to the pressure (for which Imq^ = 0), 
while the second and the third terms will be responsible 



for the EW contribution (for which Keqz = 0). It is then 
evident that the position z disappears in the Maxwell 
stress tensor. 



B. Electric and magnetic contributions to the 
pressure 

The pure electric contribution to the pressure 
Pi^'^iT, 0, 1) is due to the first term in Eq. ^ 



J^afiSaS^pv Sfj [w; ri, r2] 



ri=r2,zi=0 



, 0(1) 0(1)' 



ri=r2,zi=0 



.(1) 



El ./o \q), 'Y 



u f, \ / 



„P|2 



fe2 



(55) 



while the magnetic contribution is related by the second 
term in Eq. (|19p . and is given by 



ri=r2,zi=0 



[<9^(S33 - S22 - Sii) + Qlsii + Qls22 + QxQy{si2 + S21)] + 



[id3{QyS23 + QxSia) - id'^{QySz2 + QxSzi)]} 



D=a,zi=Z2=0 



(56) 



where s — s^^\ One can show that, as it happens for 
the equilibrium case, the magnetic contribution ([56|) co- 
incides with the electric one ([55)) . after the interchange 
of the polarization indexes s ^ p. 



C. Final expression for the pressure 



Taking the sum of (|55p and ([56]) one finds that the 
pressure P^^'^iT, 0, 1) in Eq. gS]) is 
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1 



duj . ,, ^ 



dQQ 



Req[ 



(1) 



19^ , 

+ k.p) (1 + |r^pe-«-^^') + 2 - Re (r^e^^^^')] 

^ l„(l)|2 



From this general expression one can extract the contri- 
bution of the propagating waves (PW) in the empty gap, 
for which is real and hence = kzP, and the con- 



tribution of the evanescent waves (EW), for which q^ is 
pure imaginary and hence q^ — —\qz\'^'. 



p„^eq,PW(y^ 0, I) 
pnoq.EW^y^ 0, I) 



\D. 



19^ 



Re(r^) 



ei{u;)\k 



(1)|2 



,(1)|2 



|£iM|fc2 



Re(rf) 



, (58) 
.(59) 



Now, using helpful identities [4g| and similar ones for 1 ^ 2, it is possible to express 

pncq,PW^j. Qj^ and P;;;='^'^^(T,0,0 as 

Regi'^ |tfp _ Reg4l-|rfp)+2Img, Imrf 



Re (e^(a;) gi'^) jifp ^ j^^^^ _ j^pp) ^ 2 Img, Imr'' 

|£l(^)| \qP\^ " ' 

(61) 



Note that the PW term ([62|) contains a distance inde- 
pendent contribution that will be discussed in the next 
section. 

The pressure P^^'^{0,T,l) can be obtained following 

(2) 

the same procedure but using the function s^j given by 
Eq. ([51)1 . The result can be obtained without calculation 
simply by the interchange rf^ <-^- in Eqs. (j62p and (j63p . 



I 

V. ALTERNATIVE EXPRESSION FOR THE 
PRESSURE 

The thermal pressure between two bodies in a configu- 
ration out of thermal equilibrium was derived in the pre- 
vious section, and expressed in terms of the Eqs. ((62|) and 
(I63p . In this section we present an alternative expression 
for such a pressure, explicitly in terms of the pressure 
at thermal equilibrium. In section IV Al we discuss the 
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case of bodies made of identical materials ei = S2, in 
section rVB] we discuss the general case of bodies made of 
different materials, and finally in section lV Cl we show nu- 
merical results for the pressure between different bodies 
held at different temperatures. 



A. Pressure between identical bodies 

In the case of two identical materials the pressure be- 
tween bodies can be found without any calculations us- 
ing the following simple consideration. Let the body 1 
be at temperature T and the body 2 be at T = 0, then 
the thermal pressure will be P^^'^{T, 0, 1). Because of the 
material identity the pressure will be the same if we in- 
terchange the temperatures of the bodies: P^^'^iT, 0, 1) = 
P^^'^{0,T,l). In general we know from Eq.dtl]) that the 
thermal part of the pressure is given by the sum of 
two terms each of them corresponding to a configuration 
where only one of the bodies is at non-zero temperature, 

i.e. p,";'>(ri,r2,o ^P^rm, 0,0 + nr(o,T2,0- It is 

now evident that at equilibrium, where Ti = T2 = T, the 



latter equation gives PX''(r,0,0 = P^^{T,l)/2 and we 
find for the total pressure 



(64) 



Therefore, the pressure between identical materials is ex- 
pressed only via the equilibrium pressures at Ti and T2. 
The same result was obtained for the first time by Doro- 
feyev [2^ by am explicit calculation of the pressure. It 
is interesting to note that the equation ((64)) is valid not 
only for the plane-parallel geometry, but for any couple 
of identical bodies of any shape displaced in a symmetric 
configuration with respect to a plane. 



B. Pressure between different bodies 



where we omitted the index fi. The symmetric parts, 
(1 — |rir2p) for PW and Im(rir2)/2 for EW, are re- 
sponsible for the equilibrium term P°^{T, I) /2m the non- 
equilibrium pressure as Eq. (j64p shows. Concerning the 
EW terms, if one takes the symmetric part of (|63p . one 
obtains exactly P^^'^^^'^ {T,l)/2, where Pt7^^(r,0 co- 
incides with the equilibrium EW component (jH]). The 
analysis of the PW term is more delicate, in fact if one 
takes the symmetric part (l — |7'ir2p) of (|62|). one ob- 
tains Pth ^^(T, 0/2 where 



-pCq.PW 



h 

2^ 

k 

dQQ qz 







(67) 



The above equation is different from P^^'^^ {T, I) given 
by (gll). 

The difference has a clear origin. In fact the pressure 
out of equilibrium, from which ([57]) is derived, is calcu- 
lated for bodies occupying two infinite half-spaces. On 
the contrary the equilibrium pressure Pt'h''^^(T', '^^^ 
obtained after proper regularization, and hence taking 
into account the pressure exerted on the external sur- 
faces of bodies of finite thickness [see discussion after 
Eq.dMj)]. Then the difference between (gS]) and is 
just a constant: 



-cq,PW 
th 



mo 



_ pCq.PW 



3c 



(68) 



It is convenient to present the general expression of the 
pressure in a form which reduces to Eq. (j64p in the case 
of identical bodies. It can be done using Eq. pT|) where 
P^^'^{T,OJ) is given by (EH) and and P;;^'^'^{0,T,l) 

is obtained from P^^'^{T, 0, after the interchange rf^ ^ 
'2 • 

In P^^'^{T^ 0, we can separate symmetric and anti- 
symmetric parts in respect to permutations of the bodies 
1 <-» 2. The factors sensitive to such a permutations in 
([5^ and are, respectively, 

(l-lnp) (l + |r2p) = 

(l-|rir2p) + (|r2p-|rip), (65) 

Im(ri)Re(r2) — -Im(rir2)-|- 

i [Im(ri)Re(r2) - Re(ri)Im(r2)] , (66) 



where a = n^k^/QQc^ft' is the Stefan-Boltzmann con- 
stant. Using the following multiple-reflection expansion 
of the factor \D^\^'^: 



1 \. = W I 1 + 2Re V R^e^"'"'^ | , 

|i -i?e2'9.'|2 i_ \R\2 y z^^ y ' 

(69) 

where P = r^r2 , it is not difficult to show explicitly that 
P5)) and ([^7)1 are related by The constant term in 

([68)1 comes from the first term of this expansion. 

Collecting together the symmetric and antisymmetric 
parts we can finally present the non-equilibrium pressure 
in the following useful form: 
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p: 



,eq,PW 



th 



2 

1,EW 



,oq,PW 



2 

,oq,EW 



+ 



(7^2,0 



EW/- 



This is one of the main result of this paper. Here 
B{Ti,T2) = 2a{Tf + T^) /3c is a /-independent term, 
discussed in Eq. (|68p . The equihbrium pressures 



pcq.pw^y^;^ and P^^^'^'^{T,l) are defined by Eqs.dH 



(70) 
(71) 



r 



and (I44|) and do not contain /-independent terms. The 
expressions AP^^{T,l) and AP^{T,l) are antisym- 
metric respect to the interchange of the bodies 1^2 
and are defined as: 



App^(r,/) 

AP,r(T,/) 



h 

'4^ 



duj- 



1 



,hu/kBT _ I 



dQQq, 



E 



,M|2 



,M|2 



Im [r^] Re (r^) - Im (r^) Re (r^*) 



(72) 
(73) 



Let us note that the EW term (1751) S^^s to for / 00 
because evanescent fields decay at large distances. How- 
ever, the PW term ([7^ contains a /-independent com- 
ponent since in the non-equilibrium situation there is 
momentum transfer between bodies. This /-independent 
component can be directly extracted from (|72p using the 
expansion (|69|) . This expansion shows explicitly the con- 



tributions from multiple reflections. The distance inde- 
pendent term corresponds to the first term in the ex- 
pansion (j69p . and it is related with the radiation that 
pass the cavity only once, i.e. without being reflected. 
Finally it is possible to write AP^{T, /) as the sum 
APPW(T, /) = APPW(T) + ApPJ(r, /), where the con- 
stant and the pure /-dependent terms are respectively 
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At thermal equilibrium Ti = T2 = T the sum of 
(|70| and ((7T|) provides the Lifshitz formula except 
for the term — 4crr*/3c, which is canceled due to the 
pressure exerted on the remote external surfaces of 
the bodies, as explicitly shown in the next section. 
Out of thermal equilibrium, but for identical bod- 
ies, rf^ — rj, the antisymmetric terms disappears: 
APPW(T,/) = APt|W(T,/) = 0. In this case Eq. (EH) is 
reproduced. 

It is now clear that, due to the antisymmetric terms, 
Ea. ([M)) is not valid if the two bodies are different. The 



problem of the interaction between two bodies with dif- 
ferent temperatures was previously considered by Doro- 
feyev [25|] and Dorofeyev, Fuchs and Jersch The 
authors used a different method, based on the general- 
ized Kirchhoff's law The general formalism of [2^ 
agrees with our equations ([74)) , ((75)) . However, our re- 
sults are in disagreement with the results of [2^, where 
Eq. (|M)) was found to be valid also for bodies of different 
materials. So that we argue that the results of the last 
paper were based on some inconsistent derivation. 
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C. Numerical results for the pressure between two 
different bodies out of thermal equilibrium 

In this section we show the resuhs of the calculation of 
the pressure between two different bodies, for configura- 
tions both in and out of thermal equilibrium. In figures [2] 
and [3] we show the numerical results of the pressure for a 
system made of fused silica (Si02) for the left-side body 1 
and low conductivity silicon (Si) for the right-side body 
2. In both cases the experimental values of the dielec- 
tric functions in a wide range of frequencies were taken 
from the handbook [l^]. In particular in Fig[2]we show 
the thermal pressure P^'^{Ti,T2, 1), sum of Eqs. fZO)) and 
(|7ip . as a function of the separation I between 0.5 /im 
and 5 ^m. Here we omit the ^-independent terms. The 
pressure is presented for the configuration (Ti = 300K, 
T2 = OK) [solid line] and for the configuration (Ti ~ OK, 
T2 = 300K) [dashed]. We plot also the thermal part 
of the force at thermal equilibrium, which is the sum of 
Eqs.dlSl) and (gH), at the temperature T = 300K [dot- 
ted] . The sum of the two configurations out of thermal 
equilibrium provides the force at thermal equilibrium. In 
figure [3] we show the relative contribution Pth/Po of the 
thermal component (only the Z-dependent terms) of the 
pressure with respect to the vacuum pressure Po{l) given 
by Eq.gl]). 

We performed the same analysis for a different couple 
of materials, and in particular we considered sapphire 
(AI2O3) for the left-side body 1, and fused silica (Si02) 
for the right-side body 2. Also in this case the experi- 
mental values of the dielectric functions were taken from 
the handbook fJ?] . The results of such calculations are 
shown in figures 2] and O where the same quantities of 
figure [2] and [3] were plotted. 

From figures [2] and [3] it is evident that at small sepa- 
rations the pressure at (Ti = 300K, T2 = OK) is lower 
than that at (Ti = OK, T2 = 300K), and the situation 
is inverted at large separations. This is a characteristic 
feature of the materials we use. In fact for the sapphire- 
fused silica system we found the opposite behavior, as it 
is evident from figures H] and [S] This behavior is the re- 
sult of the interplay between the relevant frequencies in 
the problem, i.e. the thermal wavelength Ay, the separa- 
tion I, and the different positions of the resonances in the 
dielectric functions for the different couples of materials. 



le-3 



T = T = 300K (dotted) 



300K,T2 = 0K (solid) 
0K,T2 = 300K (dashed) 




2 3 
separation [|lm] 



FIG. 2: Thermal component (only i-dependent part) of the 
pressure out of equilibrium for fused silica-silicon system in 
the configuration (Ti = 300K, T2 = OK) [sohd] and in the 
configuration (Ti = OK, T2 = 300K) [dashed]. We plot also 
the thermal part of the force at thermal equilibrium at T = 
300K [dotted]. 
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FIG. 3: Relative contribution of the thermal component of 
the pressure (only /-dependent part) out of equilibrium for 
fused silica-silicon system in the configuration (Ti — 300K, 
T2 = OK) [sohd], (Ti = OK, T2 = 300K) [dashed], and at 
thermal equilibrium at T = 300K [dotted]. 



VI. PRESSURE BETWEEN TWO THICK 
SLABS 

In sections IIVI and |V] we derived and discussed the 
non-equilibrium pressure between two materials filling 
two infinite half-spaces. We did not regularized the 
pressure, i.e. we did not considered the extra-pressure 
due to the presence of the external surfaces of the bodies. 
This would simply add new ^-independent terms. We 
focused mainly on the /-dependent part. In this section 
we fill this gap, and derive the exact constant terms of 



the pressure for the general case of two bodies of finite 
thicknesses at different temperatures, in presence of 
external radiation. 

At thermal equilibrium, due to the momentum's con- 
servation theorem, the pressure cannot contain constant 
terms. In fact both Eg. ([43)) and (|44)) go to zero as I 
goes to infinity. Here the regularization was performed 
by subtracting the bulk part of the full Green function 
[see discussion after Ea. ([M)) ]. The inclusion of the bulk 
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FIG. 4: Same of Fig[2j for the sapphire-fused sihca system. 
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FIG. 5: Same of FigO for the sapphire-fused sihca system. 



part would add an extra ^-independent term —AaT*/3c, 
as it is evident from the non regularized Eq. ()68|) . 
Physically the origin of this extra-term is due to the 
fact that the bodies are considered to be infinite, and 
hence, have no external surfaces. The presence of the 
external surfaces generates an extra pressure AaT'^/Sc, 
and finally the total pressure becomes Z-independent. It 
is worth noticing that at thermal equilibrium the force 
acting on one body is exactly the same (apart from the 
sign) of that acting on the second body. 

Out of thermal equilibrium, for bodies occupying 
two half-spaces, one finds the non-regularized pressure 
given by the sum of Eqs. ((^ and ([S5)) . In this case 
the pressure contains distance-independent components, 
and is the same on both materials (apart from the sign) . 
For bodies of finite thickness one should account for 
extra ^-independent terms in the pressure due to the 



presence of two more interfaces between the bodies and 
the external regions (see FiglB]) where, in general, the 
radiation is not in equilibrium with the bodies. In this 
configuration the pressure acting on the body 1 can 
be different from that acting on the body 2. It should 
be noted that the new Z— independent terms should be 
added to ([70]) and (|7T|) . and originate from the PW 
waves only. Below we derive the result for such a general 
configuration, by manipulating Eq. ()62p . 

Let us consider the case where both the bodies occupy 
thick slabs, as represented in FigEl On the left of the 
body 1 impinges radiation at temperature Tbti, while on 
the right of the body 2 impinges radiation at temperature 
Tbb2- Then the pressure acting on the body 1 and body 
2 will be respectively: 



pnoq.^ 



bbl, 7l, T2, /) 
T2, Tbb2, 



r2,0 + a(Ti,T, 



Ibi), (76) 
-P^C{Ti,T2, 1) + PKiT2,Tbb2){77) 



Here P^'^{Ti,T2,l) is the pressure out of thermal equi- 



librium given by the sum of and ([S5)) for materials 
filling infinite half-spaces. Pl is the pressure due to the 
presence of a new left-side interface of the material 1 
while Pr is the pressure due to the presence of a new 
right-side interface of the material 2. Both Pl and Pr 
are constant terms and include two contributions: the 
pressure of the external radiation impinging on the outer 
interface and the back reaction produced by the emission 
of radiation from the body to the vacuum half-space. 

For thick enough slabs, it is possible to calculate the 
terms Pl and Pr using the expression of the pressure 
acting on a body h which occupies an infinite half-space. 
In general it has a dielectric function Sh, is at temper- 
ature Tfi, and a thermal radiation with temperature Tbb 
impinges on its free surface. There are two possible con- 
figurations. One corresponds to the body h on the left 
and radiation impinging from the right, the second cor- 
respond to the body on the right and radiation imping- 
ing from the left. In the two cases the pressures can be 
expressed in terms of the pressure between two infinite 
bodies P^^'^'^^ (Ti, T2, 1) derived in the previous section, 
and are respectively: 



PR{Th,Ti 



bb) 



pr'^^^iT^^^o^i) 



£2 = 1 



PLiTh,Tbb) 



P 



ncq,PW 



tti 



{nbAi) + pr''"^io,n,i) 



£1 = 1 



(78) 



.(79) 



Here P;^^'''^'^ iT,0,l) and P;^^'''^'^ {0,T,l) are given by 
Ea. ([S^ . After explicit calculations one find 



pncq.PW 



(0,T6b,0 



£l=£h 


'^^-^ bb 


£2 = 1 


3c 


£l=£h 


2aT4 


£2 = 1 


3c 



Pd{Tbb), (80) 



Pdin), (81) 
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where 



JO Jo p 

(82) 

and are defined similar to (|35p but using the dielectric 
function e^. Finally, we obtain the main result of this 
section, i.e. Eg. ([75]) becomes 

Pn{Th,n,) = _MZ|±Im) + p^^T^) _ p,(T,,). (83) 
6c 

In the same way it is possible to calculate PhiThjTbt) 
from Eq. ((79)) . and it is evident that the result will be 



(84) 



At equilibrium — T^b — T we find that Pr(T, T) = 
—Ph{T,T) = — 4(tT^/3c does not depend on material 
characteristics and coincides with the pressure of the 
black body radiation. It is also interesting to see that 



for a white-body (W), corresponding to 



1, and 



for a black-body (B), corresponding to |r^p = 0, one 
obtains 

PR(0,r)w--^, Pr(0,T)b--^, (85) 



3c 



^r(T,0)w = 0, Pr(T,0)b 



3c 



3c 



(86) 



From these relations one can see that 
-Pr(0,T)w/-Pr(0,T)b = 2, as it should be for the 
radiation pressure. Furthermore one has that that 
Pa(T, 0)w = 0. This is the consequence of the fact that 



if 



.A* 12 



1 the radiation impinging on the surface from 



the interior of the material is fully reflected and there is 
no flux of momentum outside the body. 

In the particular case when the external radiation is at 
equilibrium with the corresponding body, i.e. Tbbi ~ Ti 
and Tbb2 — ?2, from ([55)1 and ([M]) one obtains that Eqs. 
([76)1 and (|77|) become respectively 




'bb2 



FIG. 6: Schematic figure of the two-slab system out of thermal 
equilibrium. 



VII. LONG DISTANCE BEHAVIOR OF THE 
SURFACE-SURFACE PRESSURE 

Let us consider now the surface-surface pressure in the 
limit of large separation. In this limit the relevant fre- 
quencies are uj ~ c/l <^ ksT /h. If this frequency is 
smaller than the lowest absorption resonance in the ma- 
terial, one can use the static approximation for the di- 
electrics and change ei{uj) eoi. Some dielectrics can 
have very low- lying resonances. For this case we devel- 
oped a special procedure that will be discussed later. 

At thermal equilibrium the pressure is given by 
Eas. (|43| and (|44l) for the PW and EW components, re- 
spectively. In the limit of large distances these compo- 
nents behave as 



p:^''"^{t,i) 



47r/3 ' 
AttP 



16nP 



dx 



£lO -H £20 + 1 
£lO — 1 £20 — 1 



(90) 



(91) 



Pi2iTbbi = T,,T,,T,, I) = P,r(Ti, T2, + (87) 

P2:tliTi,T2,Tbb2 = T,,l) ~P^^\T,,T2,l) - ^(88) 

If the whole system is at thermal equilibrium, Ti = T2 = 
T, the last two equations give 



P^ll{T, T, T, /) = -P^lliT, T, T, I) 



PriT,T,l) 



AaT^ 
3c 



-p:^{t,i). 



(89) 



This reproduces Eq.®, where P^^'^{T,T,l) 

■cq,PW , 



P 



th 



where C(3) « 1.2021 is the Riemann zeta-function. These 
equations are both valid at the condition 



/ > maxm=i_2 



£m0 



V £mO ^ 1 



A7 



(92) 



where At is defined in Eq. (|^ . The first term in Eq. (PT|) 
is canceled by the contribution from the propagating 
waves (|90p . and their sum provides the well known re- 
sult for the total force at equilibrium ((46|) . It is worth 
noticing that the total force at equilibrium is valid at the 
condition (|45)) . which is significantly different from ([92]) 
if one of the two bodies is rarefied. 

The surface-surface force in the non-equilibrium case is 
given by Eqs. ||7D|) and ([7T|) . Omitting the /-independent 
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terms one finds for the large distance behavior the fol- lowing result [23 |: 
I 



1 - VS20 



V £10 ~ 1 + %/ £20 



87r2Z3 



dt 



da; 



2 

X e 



1 £2oV£io — 1 + sioV £20 ~ 1_ 
ImK(t)]Re[r2^(i)] 



-2:12 ■ 



(93) 
(94) 



Here r^(t) are the Fresnel reflection coefficients (|35|) in 
the static approximation Sm = £mOj and t is defined by 
the relation = A;^ (l + t'^). Note that the equations 
([Ml and ([M]) are also valid at the condition (P^ . 

In the two following subsections IVII Al and IVII Bl we 
will describe the procedure we used to calculate the large 
distance asymptotic behaviors (l93|) and (l94|) for the PW 
and EW components, respectively. 

A. Asymptotic behavior for PW 

In this subsection we derive the expansion of the PW 
contribution P^^^"^'^"^ {T , 0, 1) at large distances, just an- 
ticipated in Eq. ((93|) . We concentrate on the ^-dependent 
part only. One can start from Eq. ((62|) . It is helpful to use 



the multiple- reflection expansion expressed by Eq. (|69p . 
The first term in this expansion corresponds to radiation 
which is emitted by one plate and absorbed by the other 
one, without being reflected back. This is a distance in- 
dependent term which we omit. All the other terms of 
the sum give contribution to the distance dependent part 
to which we are interested. 

Let us introduce new variables and parameters in 
Eq. (HI), i.e. 

^ = ^' (95) 

The limit of large distances corresponds to a <C 1. In 
terms of these new variables one finds 



^Re. 



da;- 



\^ intx / a. 



(96) 



where the reflection coefficients as functions of x and t 
are 



t - y/Tr, 



l+<^ 



t + \/Sm - 1 



t2 



Smt + V^m - 1 + ' 



(97) 
(98) 



and Sm = £m{kBTx/h) is a function of the variable x. 
For a ^ 1 the integrand in (|96p oscillates fast and it is 
possible to show that the relevant values of variables in 
the integral are a; ^ 1 and t ~ a/n. Then, expanding the 
reflection coefficients for small values of t and integrating 
over t explicitly one finds the leading term in a 



knT 



n—l 



00 



sm(na;/Q;) 

da; — 

- 1 



■E 

li=s,p 



9t.{x), (99) 



where the following functions of x were introduced 



(x) 



2Re(/?i) 
Re (/?i + P2) ' 



with 



5p {x) 



7m (a;) 



2 Re (71) 
Re (71 -1-72) 



(100) 



(101) 



The leading contribution to P^"'^'''^^(r, 0, ^) comes from 
the region a; a/n <C 1, where — \ k, x. Note that 
one can do this expansion only after explicit integration 
over t. After the change of variable y — nx/a we obtain 



r,nGq,PW 
th 



(r,o,0 



Btt^ZS E 

n—l 



.10 



E 9^.{ay/n). (102) 



y 



The relevant range of integration here is y ~ 1, and then 
the important frequencies in the dielectric functions en- 
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tering in Ea. (|10ip are of the order of ftw ~ aksT. Most 
of the dielectrics (but not all) at these frequencies have 
no dispersion in the spectrum and one can take the static 
approximation (ay/n) « (0). In this case the inte- 
gral in Eq. (|102|) can be calculated exphcitly: 



pncq.PW^y^ 0, 1) 



16nP 



C(3)[g,(0)+5p(0)], (103) 



where 



9. (0) 
5p(0) 



2V£20 - 1 



y/Sio - 1 + y/e20 - 1 ' 

2£ioV £20 — 1 

£20\/ £10 ^ 1 + EloV £20 ^ 1 



(104) 
(105) 



Eq. (jlOSp coincides with ((93|) after elementary transfor- 
mation. This expression is valid under the condition (j92[) 
that justify the expansion on i done for the reflection co- 
efficients dSll) and (|98l) . 

It is interesting to derive also the large distance behav- 
ior ((90)1 for the equilibrium case. To do this we can note 
that the symmetric part of the non-equilibrium pressure 
in respect to the interchange of the bodies coincides with 
one half of the equilibrium pressure as Eq. (j70p demon- 
strates. The symmetric part of both gs{x) and gp{x) is 
equal to 1 and we immediately reproduce the result (^0)1 . 



B. Asymptotic behavior for EW 

In this subsection we show how to evaluate the asymp- 
totic behavior of the EW contribution to the pressure 
pn^eq,EW^y^ 0, 0, whose rcsuK was anticipated in Eg. ([94)1 . 
We start from the general expression for P^^^'^"^ (T, 0, 1) 
given by Eq. (|63p . Substituting in this equation x and 

fc2 (14- ^2 V one 



a given by (|95|) . but defining < as 
finds for the pressure 

) 

dtt^e-^*/" 



E 



- 1 

Im(r^)Re(r^) 



II 



^M^M „~xt/a\2 ■ 



(106) 



Here the reflection coefficients are functions of t and x 
and take the form 



it — Ve^ 



l-i2 



+ Ve-m - 1 



i2' 



1 -t^ 



1-i^ 



(107) 
(108) 



with Em = em{kBTx/h). 

Differently from the PW component, here the relevant 
ranges of variables in the integral (|106p are a; ~ a and t ~ 
1. Small values of t do not give significant contribution 
because the integrand is suppressed by a factor t coming 
from Im (r^), that do not appears in the PW case. Then 



for large distances it is possible to expand on small values 
of X and approximate — 1 ~ x. It is convenient to 
introduce the new variable y — xt/a instead of x, for 
which the important range is now y ^ 1. In terms of y 
and t the pressure can be presented as 



ksT 



dt 







dyy 



E 



Im {r>^)Re{r!^) 



(109) 



' 1 ' 2 



The relevant frequencies in the integration are uj c/l <^ 
ksT/h and then it is possible to use the static approx- 
imation for the dielectric functions. In this approxima- 
tion the reflection coefficients depends only on one vari- 
able ''i 2(^72/) fi2W '^^^ '^^^ reproduce (after the 
change y — > x) the asymptotic behavior (j94p for the pres- 
sure Ptl^(T, 0, /) . 

The pressure in the EW sector can be presented in an 
alternative form using the multiple-reflection expansion. 
To this end one can note that 



'Xt/a 



|1 - r^r^ e-^*/"|2 Im«r^) 



Ini«r^)" -nxt/a 



(110) 



and can put this expansion in Eq. (|106p . In the static 
approximation the integral over x can be found explicitly: 



da; — 

pX 

6 



^"^*/" = «'(3)(l-hnVa), (111) 



where ^"'^^ (1 -I- nt/a) is the polygamma function [ist . 
Since a is small, one can take only the asymptotic of this 
function, which is ^''^^ (1 -I- ntja) 2{a/nt)^. Then the 
EW pressure can be presented as 



pn^cq.EW^p^ 0, 1) 



keT 



^ 1 r°° dt 



E 

lJ.=s,p 



Im(rORe(r^ )^ 
Im (r^Tj ) 



Im«r^)". (112) 



This representation is helpful for the analysis of the rar- 
efied body limit that will be presented in the next section. 

It is interesting to derive also the large distance be- 
havior (PT|) for the equilibrium case. One half of the 
equilibrium pressure, P^^" (T, l)/2, is equal to the sym- 
metric part of Eq. pi2p in respect to the bodies inter- 
change. The] 
in Eq. pi^ : 



change. Therefore, to get P^^''^'"{T, I) we have to change 



Im(r'/)Re(r^) 
Im {riv!^) 



(113) 



In this case the integrand in Eq. (|112p becomes an ana- 
lytic function of t with the poles at i = and at infinity. 
The integral can be calculated using the quarter-circle 
contour of infinite radius closing the positive real axis 
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and negative imaginary axis. This is because it Qz 
must have a positive real part. FinaUy the integral is 
reduced to the quarters of the residues in the poles and 
gives: 



,EW 



kBTC{3) 

oo ^ 



kr,T 



SttP '■ — ' n 

n—l 



£10 + 1 £20 + 1 



(114) 



The sum in this expression can be written in the equiva- 
lent integral form so that (|114p coincides with Eq. ([9T|l . 
Note that only p-polarization contributes to the pole 
at infinity. This is because at infinity r^^ — > but 
i^m (£mO - l)/(£mO + 1) stays finite. 



VIII. PRESSURE BETWEEN A SOLID AND A 
DILUTED BODY 

A case of particular interest is the interaction between 
solid and diluted bodies. In fact the first measurement of 
the non-equilibrium interaction was done between an ul- 
tracold atomic cloud and a dielectric substrate [TbS] . From 
the theoretical point of view this case is the most simple 
for analytical analysis. 

Here we investigate the pressure between a hot dielec- 
tric substrate of temperature T (body 1) and a gas cloud 
(body 2) at large distances. When the second body is 
very dilute we can consider the limit (£2 — 1) — > 0. If 
both bodies are at the same temperature T, the equi- 
librium pressure can be found by expanding Eq. (I46|) on 
small values of (£2 — 1). The leading term is 

This pressure, valid at the condition flS)) . is proportional 
to (£2 — 1) = Aima, where n is the density of the material 
2 and a is the dipole polarizability of its constituents 
(for example atoms). We can see that the pressure is 
additive since the additivity would in fact require a linear 
dependence on the gas density n and hence on (£20 — 1). 

If one performs first the diluteness limit of the exact 
surface-surface pressure, and then takes the large dis- 
tance limit, one obtains very interesting asymptotic be- 
haviors for the PW and EW contributions respectively 

M 



r,cq,PW 
th 



mo 



{ksT)^ £10 + 1 



24 /2 ch V£io - 1 
{kBTf £10 + 1 
24 P ch V£io - 1 



(£20 - 1), (116) 
(£20 - 1). (117) 



In deriving these limits we assumed that fc^T is much 
smaller than the lowest dielectric resonance of both the 
body 1 and of the atoms of the dilute body 2. Such 
asymptotic behaviors for the PW and EW components 



depend on the temperature more strongly than at equilib- 
rium and decay slower at large distances T'^/P). It is 
also remarkable that the PW component of the surface- 
rarefied body pressure (|116p depends on the dielectric 
functions and is repulsive, differently from attractive na- 
ture of the P W component of the surface-surface pressure 

The PW and EW terms pTB)l and pTZ)) exactly can- 
cel each other, and in order to find the total pressure one 
should expand the corresponding expressions to higher 
order. The final result is given by Eq. (|115p . In con- 
figurations out of thermal equilibrium there will be no 
longer such peculiar cancellations between the PW and 
EW terms. In this case the new asymptotic behavior 
~ T'^/P will characterize the total pressure at large dis- 
tances, while there will be a transition to a ~ T/l^ be- 
havior at larger distances. 

In particular, the result of the surface-rarefied body 
pressure out of equilibrium can be presented as [2^ 



keTC £10 + 1 



VI^^f{v), (118) 



where 



V£20 - 1 



(119) 



is a dimensionless variable and C = 3.83 • 10^^ is a 
constant. The function /(?;), whose expression will be 
derived below [see Eq. (ll49|) . together with Eqs. ()131|) . 
()140p . and ()14ip ]. is a dimensionless function of v. It is 
possible to show (see derivation below) that f{v) 1 for 
V 00, while f{v) — > t;/24C for v 0. This function 
is shown in Fig. [71 Eq. (jll8p is valid at the condition 
of large distances 1/Xt ^ 1 , which does not restrict the 
value of V. 

At large values of v the pressure (|118p becomes 



pnT,o,i) 



keTC £10 + 1 



1 



V£20 — 1, 



(120) 



and is proportional to -\/£2o — 1- This peculiar depen- 
dence means that the pressure acting on the atoms of 
the substrate 2 is not additive. The non additivity of 
the pressure can be physically explained as follow: for 
large I the main contribution to the force is produced 
by the grazing waves incident on the interface of the 
material 2 from the vacuum gap with small values of 
Izjk < ^/e2o ~ 1- Hence the reflection coefficients from 
the body 2 is not small even at small £20 — 1 and the 
body cannot be considered as dilute from an electrody- 
namic point of view [23| . This is a peculiarity of the non- 
equilibrium situation. In fact at equilibrium this anoma- 
lous contribution is canceled by the waves impinging the 
interface from the interior of the dielectric 2, close to the 
angle of total reflection. In a rarefied body such waves 
become grazing. Notice that the pressure (|120p is valid 
at the condition 



A7 



V£20 



(121) 
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FIG. 7: Dimensionless function f(v) [see Eqs. JTTS} and (|149) 
describing the transition between additive and nonadditive 
regimes. The dashed line presents the asymptotic limit at 
small V. 
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FIG. 8: Relevant length scales and asymptotic behaviors of 
the surface-rarefied body pressure out of thermal equilibrium. 
There is a first region given by Eq. p23p where the pressure 
is additive and coincides with Eq. (|122|l . and a second region, 
satisfying Ea. (ll2ip . where the pressure is given by Eq. p20p 
and is no longer additive. 



which becomes stronger and stronger as (^20 ~ 1) 
At small V one finds from Eq. (|118p 



(fesT)" £10 + 1 



24/2^0 Veio - 1 



(£20 - 1) 



(122) 



In this case the additivity is restored but the temperature 
dependence is not linear any more and the pressure de- 
creases more slowly with the distance. This result holds 
at distances 



At < ? < 



A7 



\/e2o - 1 



(123) 



It is worth noting that the interval (|123p practically dis- 
appears for dense dielectrics. 

The above discussion can be summarized as follow [see 
FigE] . If the dielectric 2 is very dilute but still occupies 
an infinite half space (or anyway is thick enough, in 
the sense defined above), there is a first region given by 
Ea. (|123p where the pressure is additive and coincides 
with Ea. (|122p . At larger distances, satisfying Eq. (|12ip . 
the pressure is given by Eq. (|120p and is no longer 
additive. In the intermediate region I ^ Xt/\/^20 ~ 1 
equations (|122p and (|120p are of the same order. 

It is interesting to note that, due to the diluteness con- 
dition (£20 - 1) < 1, in both regions ([^1]) and the 
thermal term APth [sum of Eqs. (|73)) and ([75)1 ] gives the 
leading contribution into the ^— dependent component of 
the total pressure P^"°'^(T, 0, 1). This clearly emerges from 
Eqs.dini) - dZI]), by comparing (fot T2 = 0) the large dis- 
tance behavior of the pressure at equilibrium P^^(T, I) 
given by Eq. ()115p . with the large distance behaviors of 
the total pressure just derived, given by Eqs. (|120p and 
(I122p . The consequences of this are remarkable. In fact 



the large distance behavior of the total pressure becomes 
proportional to 

p;r(ri, T2, ~ APth(ri, - APth(r2, 0, 

and the interaction between the two bodies will be 
attractive if Ti > T2 and repulsive in the opposite case 



Below, in sections IVIII Al and IVIII Bl we present the 
derivation of Eq. (|118p for both the PW and EW com- 
ponents, which give rise respectively to the asymptotic 
behaviors and p^ . 



A. PW contribution 

In this section we focus on the PW contribution to 
Eq. (|118p . One can do explicit calculations if the di- 
electric functions of the materials do not depend on 
frequency. This is a good approximation for the di- 
luted body 2. Since we are interested in the large dis- 
tance asymptotic, this approximation is also good for 
the solid body 1 if the material has no resonances for 
u! < c/l kBT/h. In the case of static dielectric func- 
tions the integral over x in Eq. ([TO]) can be evaluated via 
the polygamma function: 



dxx^ 



1 



pintx/a _ ^(3) ( ^ 



nt 



(124) 



Then, introducing the new variable u instead of t and the 
parameter b according to the definitions 



t 



6 = 



£10 - 1 



-\/£20 — 1 ' V £20 — 1 

one can expand ri{u) in series of 1/6 



> 1 



- 1- 



2e^u 



(125) 



(126) 
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and in the same approximation one has 



r2 



(127) 



The result is the following expression for the pressure 



r,ncq,PW 
th 



(T,0,0 = - 



2{kBTY £10 + 1 



1/Ve20-1 



l + r2 



dww3^-^(-r2)"Re*(3) (1 - i2nvu) ., 
1 ~ ''2 



(128) 



where the parameter v is given by Eq. pi9p . Here the 
sum on polarizations gave the factor £10 + 1. In the lead- 
ing approximation the integration over u can be extended 
up to infinity. Furthermore, the real part of '^^'^^ (1 — iy) 
can be presented as 



Re*(3) (1 „ ^ I A_ j^J_ _ coth^j/^ 



(129) 



After some transformations equation p28p becomes: 



ksT £10 



where fpw{v) is given by 



Ve20 - l/piv(w), 

(130) 
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FIG. 9: Function fpw{v) (solid) [Ea. (|13ip ] and its asymptotic 
limits (dashed) at small [Eq. (|135p ] and large [Eq. (|133|l ] values 
of V. 



only the n — 1 term in the sum is relevant. Then one 
obtains 



1 

fpwiy Q) ^ -—— I duux 
62tt Jq 

coth(27rvu) 



2tivu 



(134) 



and finally: 



fpw[v) = -— V — / auu 



rT^(--2)"x 

^ '2 



1 



2'Knvu 



cot]i{2Tmvu) 



(131) 



The function fpw{v) can be calculated explicitly for 
large and small values of v. When w ^ 1, the impor- 
tant range of u in the integral (|13ip is it ^ 1 and one 
can expand the reflection coefficient r2 on small values of 
u. Then the function fpw{v) is reduced to 



00 



fpw{v 00) = 2^ — 



167r 



1 



duu X 



2'Knvu 



coth(27rnuw) 



(132) 



The integral here is easily calculated by parts and finally 
one finds 



fpw{v 



Cpw ~ 



C(3) 



(133) 



!pw{v ^ 0) = — . 



(135) 



The function fpwW) ^^nd its asymptotic behaviours at 
large and small v are shown in Figl9l 



B. EW contribution 

The derivation of the EW component of the pres- 
sure (|118l) can be performed starting from the expres- 
sion (|106p for the pressure 

p^„oq,EW(j.^Q^^) By perform- 
ing the multiple-reflection expansion with the help of 
Eq. (|110p and calculating the integral over the variable x 
using piip one finds that Eq. (|106p becomes 



p-"-(T,o,o=^f:rd«^x 



_i ^0 
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^ Im(r5'r^) 
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(136) 



When t; ^ 1 the significant values of u in the integral 
(I13ip are w 3> 1, and one can make the corresponding 
expansion in the reflection coefficient (|127p . In this case 



As in the case of the propagating waves one can introduce 
the variable u instead of t according to Eq. (|125p , and can 
make the expansion for large h. Then for the reflection 
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coefficients one gets 
^ i2u 



i2siu 



(137) 



Now one should distinct the integration ranges < u < 1 
and 1 < u < cxD, since the integrands are different in 
these ranges. Let us do it with the superscript (1) or (2), 
respectively. 

As in the case of propagating waves (|130p the pressure 
can be presented as a parameter-dependent factor times 
a universal function of v ~ 1^/s2q-^/^t' 



PX"*' (r, 0, /) 



(138) 

where the function fEw{v) includes contributions from 
< M < 1 and 1 < u < oo ranges: 



(139) 



For these functions one has the following expressions 



lEW 
2^2-1 



n=l 
j3 







, ^Im(-r2)"-3-^*(l + 2nTO); (140) 



1 1 /'OO 

n=l ^ 



{-r2r X 



— *(l + 2nw). (141) 
dif^ 

Here we have used the relation between the polygamma 
functions 



(142) 



where (1 + y) is the digamma function [48*1. 

Let us discuss now the asymptotic behavior of the func- 
tions fEw'^\'^) small and large values of v. For large 
V the contribution from the range u < 1/w in the integral 
Eq. (|140[) is negligible, and one can consider the digamma 
function at large arguments ^'(1 -f 2nvu) ln(2nuu). 
Then the integral can be calculated after the substitu- 
tion u = siniy9. It gives the following result 



1 1 

n=l 



+ (_l)W^(„ + 3/2) + 



where the /3-function is defined as 



4^2 - 1 



l + y 



*(l) 



(143) 



(144) 



(2) 

To find f^wi''^) '"^^ ^^^^ asymptotic value 
of \l/(l -I- 2nvu), make the change u — coshx, and after 
the integration one obtains 



—y 



2("1)^ 



47r2 



^ n^{^.■n?■ - 1) ' 



(145) 



Taking the sum of both functions /^j^^ and f^^^ one finds 
finally the large v asymptotic for jEwiy): 

1 °° 1 

n=l 

I + (-!)" (^(" + 3/2) + !^)]. (146) 

This sum is just a number equal to 

fEwiv oo) = Cew = -0.96 • 10-2. (147) 

Combining together the large v contributions from PW 
(fT33| and EW (fT47| one can find the constant in 
Eq. miD, i.e. C = Cpw + Cew = 3.83 • IQ-^. 

In the limit of small v it is not difficult to show that 
fswi'^) ~ and can be neglected. The main contri- 

(2) 

bution to fEw^"^) comes from the range u ^ 1/v ^ 1. 
For these values the refiection coefficient r2(u) ~ 1/4^2 
is small and only the n = 1 term in the sum is relevant. 
Then the integral over u can be calculated by parts and 
one obtains 



fEwiv 0) 



V 

48' 
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Let us note that for small the PW and EW contribu- 
tions coincide. 

The function fEwi^) shown in Fig[TOl The inset 
demonstrates the cubic behavior at small v. It should 
be noted that /ew^''^) fpwiv) approach the large v 
asymptotics rather slowly, but the sum of these functions 
reaches the large v limit faster as FiglT] demonstrates. 
The function /^^^/(f) is presented in Fig[TlJ One can see 
that it behaves in accordance with expected asymptotics. 

Finally, one can establish the correspondence between 
the function f{v) entering the general formula pi8p for 
the pressure in the limit of one diluted body and the func- 
tions fpw{v), fswi'")' ^^'^ fswi'") given by Eqs. p3T|l . 
(|140p . and (|14ip . respectively. This correspondence is 
given by the simple relation 

C fiv) - fpwiv) + fl^liv) + .f%{v). (149) 



IX. LARGE DISTANCE BEHAVIOR OF THE 
SURFACE-ATOM FORCE OUT OF THERMAL 
EQUILIBRIUM 

It is interesting to recover the asymptotic results of the 
surface-atom force out of thermal equilibrium (obtained 
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FIG. 10: Function (solid) [Eq.dUO}]. The asymptotic 
limit at large values of v [Ea. (|143|l ] is shown by the dashed 
line. The inset demonstrates v'^ behavior at small v. 




FIG. 11: Function fswiv) (solid) [Ea. lfTiTjl ] and its asymp- 
totic limits (dashed) at small and large [Ea. p45|l ] values of 

V. 



in |17| ) from the general expression of the pressure given 
by Eqs. ([S^ and To do this it is crucial to carry out 
the limit (£2 ~ 1) = 47rna2 before taking the limit 
of large distances. To show this, let us focus first on 
the EW term given by Eq. ([63]) , and perform the rarefied 
body expansion (body 2) assuming that ^£20 — 1 is the 
smallest quantity, also with respect to \qz\/k. Due to the 
effect of the Bose factor, only the frequencies u ^ ksT / h 
are relevant in the integration, and due to the exponential 
g-2(|9z| ^jjg j-gigvant wave- vectors are given by 



\qz\/k - At A > \/e20 - 1- 



(150) 



In this way at large distance it is easy to reproduce the 
Eqs.(10)-(11) of [l7|: 



^(£20 - 1) 

P 87r2 c 



dw 



\e,{uo) - 1| + [eU^) - 1] 



ei{uo)~l 



II 



• (151) 



%/2|ei(L^) 

In deriving Eq. (|15ip we also replaced £2 (w) with its static 
value £20, which is reasonable if ksT is much smaller than 
the lowest atomic resonances, and also ensures that the 
atoms of the dilute body 2 cannot adsorb the thermal 
radiation. 

For a rarefied body one has that £20 — 1 ~ 47raona, 
where Ua is the number of atoms of the body per unit 
volume and is the static polarizability of an atom. 
The pressure in this case is proportional to Ua and the 
force acting on an individual atom can be calculated as 



[^neq.EW 

th — 



n„ dl 
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It is easy to check that substituting Eg. pSip into (|152p 
one obtains exactly Eq.(lO) of Tz*]. 

However, there is also the PW contribution. The ex- 
pansion in the I -dependent part of the PW pressure 
([S^ produces a contribution identical to the EW one, 
thereby doubling the value of the force (|152[) . This ap- 
parent contradiction can be easily solved by the following 
arguments. The problem approached in the present pa- 
per is not equivalent from that approached in Ref. [17| . 
Here we assume that the second slab, being rarefied, is 
still thick enough to absorb black body radiation from 
the first slab. On the contrary, the transition to indi- 
vidual atoms (which is the case discussed in Ref. flT\) 
demands to completely neglect the absorption. Then, to 
calculate the surface-atom force correctly, one must con- 
sider the limit £3 ^ at finite thickness L of the slab 2. 
On the contrary using the expression (|62p means taking 
the opposite limit procedure, i.e. first L — s- 00 and later 
£2 0. The reason why the first limiting procedure is 
correct in this case, is that if the slab 2 does not absorb 
radiation completely, one should also take into account 
the pressure acting on the remote surface (i.e. the ex- 
ternal one), generated by the radiation coming from the 
left. In absence of absorption it is possible to show that 
the inclusion of the remote surface in the slab 2 results 
in a relatively small value of the PW pressure. Details of 
calculations are presented in the Appendix [Bl We only 
notice here that neglecting of absorption actually requires 
the condition £3 <C X^/IL. 

As a consequence, for a finite slab of rarefied gas with- 
out absorption the EW contribution (|15ip provides the 
total pressure and is equivalent to equations (lO)-(ll) of 
(17i] for the surface-atom force. In particular at tempera- 
tures less than the lowest resonance in £i(w) the pressure 
(|15ip (and hence the total pressure) takes the form [l3| 



(kBTf £10 + 1 



(£20 - 1). (153) 
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The above result holds at distances (|123p . 



X. CONCLUSIONS 

In this paper we generalized the Casimir-Lifshitz the- 
ory for the surface-surface pressure to a situation out of 
thermal equilibrium, when two bodies are kept at differ- 
ent temperatures in a stationary configuration. In con- 
trast with the equilibrium case, the non-equilibrium force 
cannot be presented as the sum over imaginary frequen- 
cies and one has to work in the real frequency domain. 
At real frequencies it is natural to separate contributions 
from propagating and evanescent waves. The delicate 
interplay between these contributions set the total force. 

For bodies made of similar materials the pressure is 
expressed via the forces at equilibrium. In the general 
case there is an additional contribution to the pressure, 
which is antisymmetric in respect to interchange of the 
materials. The propagating part of the force contains 
distance independent terms, due to the presence of an 
energy flux between the bodies in absence of equilibrium. 

We presented a detailed analysis of the force, with par- 
ticular attention for large separations/high temperature 
behaviors. At equilibrium significant cancellations be- 
tween PW and EW contributions occur. Such cancella- 
tions are less pronounced in the non-equilibrium situa- 
tion. It is established that at large distances the force 
between heated (T) and cold {T = 0) bodies behaves 
similar to the Lifshitz limit, ~ T/l^, but with different 
numerical coefhcient. However, this result is true only 
for dense bodies. If one of them is diluted the behavior 
of the force can change. 

Special attention was devoted to the case when one 
body is diluted. This is an important situation from 
which one can recover the interaction between a body 
and a single atom. Two remarkable results are found 
for this situation [23[. First, at very large distances, 
I ^ At/-\/£20 — 1, the pressure becomes non-additive, 
in contrast with the equilibrium case. Namely, the non- 
equilibrium pressure is proportional to the square root 
of the density of the diluted body, while in the equilib- 
rium it is proportional to the first power of the density 
and, therefore, it is additive. The second result concerns 
smaller distances, Xt I Xt/V^w — 1- In this case 
we found a new asymptotic behavior for the pressure, 

T'^/P, that decays with the distance more slowly than 
the Lifshitz limit at equilibrium, and has a stronger tem- 
perature dependence. A careful analysis of the transition 



region between these two limits was done both analyti- 
cally and numerically. 

The pressure between diluted and dense bodies in the 
distance range At ^ ^ <C Xt/\/£20 — 1 is used to deduce 
the surface-atom force. Earlier and with different meth- 
ods it was found in that at large distances this force 
must behave as ~ T^/Z^. The direct transition from the 
case of the surface-diluted body provides a force which 
is two times larger than that in Ref. [l3|, and both EW 
and PW terms contribute in the same way. We provided 
a detailed explanation why if the atom does not absorb 
radiation one has to neglect the contribution of the PW 
term, hence recovering the known result. 
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APPENDIX A: GREEN FUNCTIONS FOR TWO 
PARALLEL DIELECTRIC HALF-SPACES 

In this section we present the Green function, which is 
a solution of Eq.([7]). We use the Sipe Green- function for- 
malism [37| for surface optics. Sipe formulated the prob- 
lem in terms of s— and p— polarized EM vectors waves, 
in of the Fresnel coefficients of the interfaces. Here we 
use the lateral Fourier transform representation for the 
Green's function: 



(Al) 

In our geometry the Fourier transform gij [w; Q, z, z'] de- 
pends only from the modulus Q — |Q|. 



1. Green's function with the source and the 
observation points in the vacuum gap 

If both the observation point r and the source point r' 
are in the vacuum gap, the Green function can be written 
as the sum Gy [w; r, r'] = G^^=- [w; r, r'] + Gj'/- [a;;r,r'], 
of a scattered and bulk part. In particular the Fourier 
transform of these terms are 1391] : 
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gl]- [w; Q, z, z'] 



2mk^ ^ 1 



(A2) 



5^"- [c.;Q,z,z'] 



-47r(5i3<5j3'5(z - z') + 
27riA:2 



?z '- 



(z~z') 



[z - z') + e^,,(-) e^,,(-) e-<^^(--') - z) . (A3) 



Here the multiple reflections enter only in the scattered the source point r' in the body 1 or 2, are respectively 
term and are described by the denominator [s^ : 

= 1 -r^'r^e^'"-'. (A4) 



2. Green's function with the source in a body and 
the observation point in the vacuum gap 

The Fourier transform of the transmitted Green func- 
tions with the observation point r in the vacuum gap and 



-,(1) 



Hz fi—s,p 



(2) 

Qz fl—S.p 



(A5) 
(A6) 



The symmetry of the problem becomes clear when one 
set the origin of the coordinate axis in the center of the 
vacuum gap, by changing z — > z — Z/2 and z' z' — Ijl 
in Eq.j^, and ([MI- 



APPENDIX B: FORCE ACTING ON A 
RAREFIED SLAB 

As discussed in Sec. IIXI in order to recover the surface- 
atom force starting from the surface-surface expression, 
one must consider the rarefied body as occupying a slab 
of finite thickness. In this case, for non absorbing atom 
the PW term of the pressure is negligible, and the EW 
one reproduces entirely the surface-atom force derived in 
[TtI . In this section we discuss this problem and show 
explicitly that the PW term can be neglected. Let us 
consider the problem of the thermal forces between a 
body 1 at temperature T, which occupies the half-space 
(z < 0), and a body 2 at zero temperature which occupies 
a slab of thickness L in the region (I < z <l^L). In the 



gap < z < Z (region 0) and outside of the slab z > l^-L 
(region 3) we can take e = 1. The force per unit of area, 
acting on the slab in z— direction, is 

P (T, 0) = P(") - P(3) ^ ^y(o) ^ _ ^y(3) Y (Bi) 

where Tf^ and Tf) are the zz— component of the 
Maxwell stress tensor in vacuum, calculated in the re- 
gions and 3, respectively. For a completely absorbing 
slab there is no field in the region 3, Tzz — and one 
return to Eq.(IT]). Of course from (jBl[) one can calculate 
the force acting on a slab of arbitrary thickness, and can 
recover the results of this paper relative to a thick slab. 
Here we we assume that the slab is rarefied: 

£20 - 1 « 1 • (B2) 

Our goal will be to prove that for a slab without ab- 
sorption the propagating waves give the contribution 
pPW ^ pEW^ hence can be neglected. For the proof 
it is enough to consider a monochromatic component of 
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the thermal radiation impinging on the surface of the 
body 2 with the wave vector k and polarization ^ — s,p. 
In the terms of the complex amplitudes of the fields its 
contribution to the pressure can be written as [we omit 
(/i, k) arguments of the fields] 



(B3) 

The fields in the region are the sums of incident (+) 
and reflected (— ) waves: 

E(") = e(o+' + e(o-), H(o) = H("+) + H("-) , (B4) 

where H(o+) oc e*«-^ and E(o-),H(o-) oc e-*'?-^ 

An important point of the proof is that incident and re- 
flected waves give independent contributions to the stress 
tensor: 



7^(0) ^ 7^(0+) I rp{o- 



(B5) 



The additivity property (|B5p is obvious. Presence of the 
mixed term containing both E*^°+^ and E^°~)* would re- 
sult in the z— dependence of Tzz- But this is not possible 
since it violates the momentum conservation. 
By definition we have that 



E(°-) k) = E(°+) (/i, k) 



(B6) 



where _r(^''') is the reflection coefficient from the slab, 
for the (/i, k)-wave . Taking into account the Fresnel 
relations between the field components at the reflection, 
we easily find that 



rir) = i?ri°+), Ti? = (i + i?)Ti°+). 



(B7) 



Let us consider now the fields in the vacuum region 3. 
There is only a refracted wave and we have |E*^'^)| = 
£)ip-M) 1^, where ZJ^'*''^) is the transmission coeffi- 

cient. In absence of absorption D'^'^'^^ = 1 — R'-^^'^K This 
means that 



Ti3)(^,k) = (l-i?(^'>^))Tir)(^,k) 

1 - 



1 + i?(A',k) 



and from (IBll) one has 



P^^(M,k) 



2i?(M,k) 



1 + 



(B8) 



(B9) 



One can easily calculate R'^f^'^\ (see, for example, the 
problem N.4 in § 66 of [30|]). At real £20 ^ 1 one gets, 
independently on the polarization, the result: 

, . . sin^ T— cos So] , ,2 , ^ 

^1"'^'- ' . .(^^0) 

where is the angle of incidence. This equation is valid 
at the condition cos 6*0 ^ — 1 • Let us note that 

the surface-atom force equations of [17] must be valid 
in the "additive" regime of the Sec. IVIIIi where just 
the incident angles cos 6*0 = Qz/k ^ Xt/1 ^ V^20 — 1 
are important [see Eq. (|150p ]. For such angles R^^^'^'> ^ 
{We2Q - 1/At)^ < 1, and P^w ^ ^^^^ 
assumed that I Xt- For I < At one gets sim- 
ply P(''''^) - (£20 - 1)^- It is not difficult to check 
that ri? - (£20 - 1) P^^. Finally we find that 
pPW ^ pEW g^^^ hence the propagating waves contri- 
bution can be neglected. 

Let us discuss now the role of a weak absorption. Con- 
sider the case 



£2=£2+i£2,£2 <l,e2~l• 



(Bll) 



It is not difficult to generalize (jB9[) for a slab with ab- 
sorption: 



D 



1 + R 



rpiO) 



(B12) 



where the transmission coefficient D < 1 — R. IfP<Cl, 
pPW« (l-Pi)ri?. (B13) 



According to Problem 4 of § 66 in |3(>], one has that 
D '-^ exp ^— ^"gg ^ . The imaginary part £2 (w) must be 
taken in this estimate at uj ksT /h. The factor (1 — D) 
and correspondingly pPW small if u)Le'2 <C ccos^o- 
This gives the condition for neglecting the absorption: 



£2 (w ~ ksTlh) < \1/IL. 



(B14) 



Let us note that for evanescent waves Piz"*^^ does not 
depend on z, while the field of an evanescent wave goes 
to zero at z — !■ 00. This means that TzV^^ = 0. 
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